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ABSTRACT 

In this series of papers we investigate the orbital structure of 3D models representing 
barred galaxies. In the present introductory paper we use a fiducial case to describe all 
families of periodic orbits that may play a role in the morphology of three-dimensional 
bars. We show that, in a 3D bar, the backbone of the orbital structure is not just the xl 
family, as in 2D models, but a tree of 2D and 3D families bifurcating from xl. Besides 
the main tree we have also found another group of families of lesser importance around 
the radial 3:1 resonance. The families of this group bifurcate from xl and influence 
the dynamics of the system only locally. We also find that 3D orbits elongated along 
the bar minor axis can be formed by bifurcations of the planar x2 family. They can 
support 3D bar-like structures along the minor axis of the main bar. Banana-like orbits 
around the stable Lagrangian points build a forest of 2D and 3D families as well. The 
importance of the 3D xl-tree families at the outer parts of the bar depends critically 
on whether they are introduced in the system as bifurcations in z or in i. 

Key words: Galaxies: evolution - kinematics and dynamics - structure 



1 INTRODUCTION 



X 



A thorough understanding of the orbital structure in a 
barred galaxy potential can provide useful insight to the stel- 
lar dynamics of barred galaxies and therefore to the dynam- 
ical evolution of these objects, as reviewed e.g. by Athanas- 
soula (1984), Contopoulos & Grosb0l (1989), Sellwood & 
Wilkinson (1993) and Pfenniger (1996). Stable periodic or- 
bits trap around them regular orbits and thus constitute 
the backbone of galaxy structure (Athanassoula, Bienayme, 
Martinet et al. 1983). Thus the appearance of a given mor- 
phological feature can often be associated with the proper- 
ties of one of the main families of periodic orbits. In the 
'90s, starting with Athanassoula (1992a; 1992b), many pa- 
pers have pointed out that the gaseous response to steady 
barred potentials is, to a large degree, determined by the 
morphology of the periodic orbits in the corresponding stel- 
lar case. Thus, orbital and gaseous dynamics are linked. This 
has provided added incentive for studies of the morphology 
and the stability of periodic orbits in Hamiltonian systems 
representing disc galaxies. 

Orbital theory has often provided useful information on 
the structure of galactic bars. Thus it is now understood that 
a bar is basically due to regular orbits trapped around the 
so called 'xl' pe riodic orbits, which are elonga ted along the 
bar major axis (Contopoulos & Grosb0l 1989). Such orbits 
do not extend beyond the corotation resonance, and in many 



cases no suitable elongated orbits can be found beyond the 
4:1 resonance. This led orbital the ory to predict that bars 
should end at or before corotation (Contopoulos 1980). Or- 



bital theory was also able to predict - at the right distance 
from the center - the loops of the near-infrared isophotes (see 
the case of NGC 4314 in Patsis, Athanassoula & Quillen 
1997). Yet not all important morphological features have 
been so far explained with the help of periodic orbits. Thus 
orbital theory has difficulties to explain the boxy isophotes 
surrounding the bars of, mainly, early-type barred galaxies 
(Athanassoula 1996, Patsis et al. 1997). Another point still 
under discussion is the morphology of the peanut-shaped 
bulges observed in edge-on disc galaxies. They are consid- 
ered by many authors as revealing the presence of a bar, 
and to be associated with the 2:1 vertical resonance. It is 
not clear, however, which families can make this vertical 
structure. Could a bar without a vertical 2:1 resonance be 
boxy or peanut-shaped when viewed edge-on? Could we have 
stellar rings out of the equatorial plane at the ILR region? 
Furthermore, the detailed dynamics of the corotation region 
and the differences in the vertical structure between fast and 
slow bars remain open issues. 

In this series of papers we use orbital theory to address 
the above questions. This is a first step towards understand- 
ing both the orbital behavior in A^-body models and the 
responses of gaseous discs to potentials derived from near- 
infrared observations. The differences between our model 
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and the well studied corr esponding 2D case of the Ferrers 
bar ( Athanassoula 1992a) reflect the changes due to the in- 



clusion of the third dimension. In separate papers we address 
the question of the morphology of the peanut shaped-bulges 
(Patsis, Skokos & Athanassoula, 2002a - paper III) and of 
the boxy isophotes of bars seen face-on (Patsis, Skokos & 
Athanassoula, 2002b - paper IV). 

Our first goal is to make a thorough study of the or- 
bital structure in 3D barred potentials, to classify the im- 
portant families, and to follow their morphological evolution 
as a function of the Jacobi integral. We start with a fiducial 
case. Many of the families we find in this model have been 
previously mentioned (e.g Heisler, Merritt & Schwarzschild 
1982; Pfenniger 1984, 1985b; Martinet & de Zeeuw 1988; 
Hasan, Pfenniger & Norman 1993). However, other, equally 
important families, have not yet been studied. 

Studying the orbital stability in a Hamiltonian system 
approximating the dynamics of a barred galaxy we get the 
periodic orbits that could be used as building blocks for a 
density model. The general rule is to look for stable periodic 
orbits, since they trap around them the regular orbits. Not 
all of them, however, are equally important. The isodensities 
of the model we use show us the topological limits within 
which we should look for the significant orbits. Stable repre- 
sentatives of families of periodic orbits which do not support 
the imposed morphology, i.e. that of a bar embedded in an 
axisymmetric disc with a central bulge, should be consid- 
ered as less important. As an example, in 2D models, let 
us mention the case of the retrograde family x4, which is 
stable over a very large interval of energies (Athanassoula 
et al. 1983). This family, however, should get a minimum 
weight when one tries to construct a self consistent model 
using a Schwarzschild (1979) or a Contopoulos & Grosb0l 
(1988) method. In a model of a 3D disc galaxy, besides the 
counter-rotating x4 family on the equatorial plane, one has 
to filter out also stable orbits with large \z\, i.e. orbits with 
large mean vertical deviations, since these orbits do not con- 
tribute much to the density of the barred galaxy. 

This paper is organized as follows: In section 2 we re- 
view briefly the parts of orbital theory that are necessary 
for understanding this paper. In particular we explain the 
use of characteristic and stability diagrams in following the 
dynamical evolution of a family of periodic orbits. We de- 
scribe also the various types of instabilities encountered in 
3D Hamiltonian systems and we introduce the nomenclature 
of the main families. The latter is necessary since a number 
of the families presented here have not been previously dis- 
cussed and thus need to be incorporated in a unique nomen- 
clature scheme. In section 3 we introduce our 3D model and 
the orbital structure in a 2D counterpart. In section 4 we 
present the main families xl, x2 and x3 and their bifurca- 
tions. In section 5 we describe the orbits around L4 (or L5) 
and around Li (or L2), as well as families outside corotation. 
We conclude in section 6. 



2 A SHORT INTRODUCTION TO PERIODIC 
ORBITS IN THE PRESENT CONTEXT 

2.1 Periodic orbits and their stability 

In this section we will briefly review some parts of orbital 
theory which are necessary for the understanding of this 



paper. A clear, easily readable introduction to the subject 
has been given by Sellwood & Wilkinson (1993). We also 
refer the reader to the pioneering works of Pfenniger (1984, 
1985b) and Contopoulos & Magnenat (1985). 

We study the stability of simple-periodic orbits in a 
barred potential in cartesian coordinates. The 3D bar is ro- 
tating around its short z axis. The x axis is the intermediate 
and the y axis the long one. The system is rotating with an 
angular speed fit and the Hamiltonian governing the motion 
of a test-particle can be written in the form: 

H = -[pi +P 2 y +pl) + V(x,y,z) - il b (xp y - yp x ), (1) 

where p x , p y , and p z are the canonically conjugate momenta. 
We will hereafter denote the numerical value of the Hamil- 
tonian by Ej and refer to it as the Jacobi constant or, more 
loosely, as the 'energy'. The corresponding equations of mo- 
tion are: 



x =p x + D. b y, y = p y 
dV n 

Px = --^+flbPy, Py 
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The space of section in the case of a 3D system is 4D. 
The equations of motion are solved for a given value of the 
Hamiltonian, starting with initial conditions (xo,xo, zo, zo) 
in the plane y=0, for y > 0. The next intersection with 
the y=0 plane with y > is found and the exact initial 
conditions for the periodic orbit are calculated using a New- 
ton iterative method. A periodic orbit is found when the 
initial and final coordinates coincide with an accuracy at 
least 10~ 10 . The integration scheme used was a fourth order 
Runge-Kutta scheme. 

The estimation of the linear stability of a periodic or- 
bit is based on the theory of variational equations. We first 
consider small deviations from its initial conditions, and 
then integrate the orbit again to the next upward inter- 



section. In this way a transformation T : R 



R 



established, which relates the initial with the final point. 
The relation of the final deviations of this neighboring orbit 
from the periodic one, with the initially introduced devi- 
ations can be written in vector form as: £ = M£o- Here 
£ is the final deviation, £0 is the initial deviation and M 
is a 4 x 4 matrix, called the monodromy matrix. It can be 
shown that the characteristic equation is written in the form 
A 4 + a\ 3 + f3X 2 + a\+ 1 = 0. Its solutions = 1, 2,3,4) 
obey the relations Ai A2 = 1 and A3 A4 = 1 and for each pair 
we can write: 



Ai,l/A< = 3 [-&i±(&?-4)i], 



(3) 



where h = 1/2 (a ± A 1/2 ) and A = a 2 - 4Q3 - 2). 

The quantities 61 and 62 are called the stability indices. 
If A > 0, |&i| < 2 and I&2I < 2, the 4 eigenvalues are on 
the unit circle and the periodic orbit is called 'stable'. If 
A > 0, and |&i| > 2, |6 2 | < 2, or \b 2 \ > 2, |&i| < 2, two 
eigenvalues are on the real axis and two on the unit circle, 
and the periodic orbit is called 'simple unstable'. If A > 0, 
1 61 1 > 2, and &2 > 2, all four eigenvalues are on the real axis, 
and the periodic orbit is called 'double unstable'. Finally, 
A < means that all four eigenvalues are complex numbers 
but off the unit circle. The orbit is characterized then as 
'complex unstable' (Contopoulos & Magnenat 1985, Heggie 
1985, Pfenniger 1985a,b). We use the symbols S, U, D, A 
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to refer to stable, simple unstable, double unstable and 
complex unstable periodic orbits respectively. For a general 
discussion of the kinds of instability encountered in Hamil- 
tonian systems of N degrees of freedom the reader may refer 
to Skokos (2001). 

The method described above has been initially pre- 
sented by Broucke (1969) and Hadjidemetriou (1975), and 
has been used in studies of the stability of periodic orbits in 
systems of three degrees of freedom. The reader is referred 
to Pfenniger (1984) and Contopoulos & Magnenat (1985) 
for an extended description. 

A diagram that describes the stability of a family of pe- 
riodic orbits in a given potential when one of the parameters 
of the system varies (e.g. the numerical value of the Hamil- 
tonian Ej), while all other parameters remain constant, is 
called a 'stability diagram' (Contopoulos & Barbanis 1985, 
Contopoulos & Magnenat 1985). With the help of such a 
diagram one is able to follow the evolution of the stabil- 
ity indices b\ and 62, and the transitions from stability to 
instability or from one to another kind of instability. We 
will loosely refer to the b — 2 and 6 = —2 lines on a sta- 
bility diagram as the b — 2 and b = — 2 axes. The S— >U 
transitions, when one of the stability indices has an inter- 
section with the b — —2 axis, or tangencies of the stability 
curves with the b — —2 axis, are of special importance for 
the dynamics of a system. In this case a new stable family 
is generated by bifurcation of the initial one and has the 
same multiplicity as the parent family. That means that the 
periodic orbits of the bifurcating family have, before clos- 
ing, as many intersections with the plane y—0, for y > 0, as 
the orbits of the parent family. The new family may play an 
important role in the dynamics of the system. S — >U transi- 
tions after an intersection of a stability curve with the 6 = 2 
axis, or tangencies of a stability curve with the b — 2 axis, 
also generate a stable family but are accompanied by period 
doubling. This means that the number of intersections with 
the plane y=0 (always with y > 0), needed for the periodic 
orbits to close, is double the corresponding number of the 
parent family. Since the most important families we examine 
here are simple-periodic, i.e. of multiplicity 1, intersections 
or tangencies of their stability indices with the b — 2 axis 
introduce in the system families of orbits with multiplicity 2. 
U— >D and D— >A transitions do not bring new stable families 
in the system and thus in principle are only of theoretical 
interest. As we will see, however, the evolution of a family 
which is found to be initially unstable may be very impor- 
tant for the dynamics of our model. The family could simply 
become stable in another energy interval, or it may play a 
major role in a collision of bifurcation s, an inverse bifur - 
cation or other dynamical phenomena (Contopoulos 1986). 
Finally in the case S— >A we have in general no bifurcating 
families of periodic orbits. 

An other very useful diagram is the 'c haracteristic' di- 
agram (Contopoulos & Mertzanides 1977). It gives the x 



coordinate of the initial conditions of the periodic orbits of 
a family as a function of their Jacobi constant Ej. In the 
case of orbits lying on the equatorial plane and starting per- 
pendicular to the x axis, we need only one initial condition, 
x, in order to specify a periodic orbit on the characteris- 
tic diagram. Thus, for such orbits this diagram gives the 
complete information about the interrelations of the initial 
conditions in a tree of families of periodic orbits and their bi- 



furcations. However, even for orbits completely on the equa- 
torial plane, but not starting perpendicular to the x axis 
we need to give initial conditions as position-velocity pairs 
(x, x) and the characteristic diagram is three-dimensional 
(Ej,x,x). In the general case of orbits in a 3D system, one 
has a set of four initial conditions and the characteristic 
diagram is five-dimensional. The representation of such a 
diagram is difficult, but when necessary we will give just the 
(Ej,x) projection. (Ej,x) diagrams that can be compared 
with the corresponding 2D models will always be given. In 
all characteristic diagrams the region to which the orbits are 
confined is bounded by a curve known as the zero velocity 
curve (ZVC), since the velocity on it becomes zero. 

2.2 The nomenclature of the main families 

Our orbital study is more extended than previous ones 
and thus brings in new families of orbits which have not 
been studied so far. We were thus brought to introduce a 
new nomenclature system, extension of the Contopoulos & 
Grosb0l (1989) system, which covers all the new types of 
orbits. 

For the main 2D families of simple periodic orbits the 
nomenclature in the present paper follows the standard no- 
tation of Contopoulos & Grosb0l (1989). We thus have the 
xl family, where orbits are elongated along the bar and 
which is the main family in the case of barred potentials, 
families x2 and x3, whose orbits are elongated perpendic- 
ular to the bar, and the retrograde family x4. 2D families 
bifurcated from xl at the 3:1 resonance region on the equa- 
torial plane are denoted by tl, t2, for consistency with 
the names used in Patsis et al. (1997). 2D families bifurcated 
at the 4:1 resonance region on the equatorial plane are called 
ql, q2, q3,.... Planar orbits related with the 1:1 radial res- 
onance will be called ol, o2.... They are encountered only 
in some models. The fiducial case presented in the present 
paper is not one of them. 

Further planar families appear beyond the xl family, at 
the gaps of the even resonances 4:1, 6:1, 8:1 etc. They are 
given the names 'f, 's', 'e'... respectively. These families, not 
directly related to the morphological problems we address 
in this series, will be discussed elsewhere. 

We name the 3D families bifurcated from the basic fam- 
ily xl at the vertical resonances as xlvn, where n denotes 
the order of their appearance in our fiducial model A (see 
below section 4). This is a convenient model to be used for 
our nomenclature, since there are families of 3D orbits as- 
sociated with all basic vertical resonances. So xlvl is the 
one bifurcated at the first S— >U transition, which happens 
at the vertical 2:1 resonance region, xlv2 is the one bifur- 
cated at the U— +S transition (second stability transition of 
the model also at the vertical 2:1 resonance region), xlv3 
is the one bifurcated at the S— >U transition at the verti- 
cal 3:1 resonance and so on. Further bifurcations of these 
xlvn families are indicated with an l .n' (for the n-th bifur- 
cation) attached to the name of the parent family; i.e. the 
first bifurcation of xlvl will be xlvl.l, the second xlvl. 2, 
etc. Further bifurcations of these families will be indicated 
by further '.n' attached to the name of the parent family. 
Thus xlvl. 1.1 is the first bifurcation of xlvl.l. The naming 
system is thus extendable at will. 

In general at each vertical resonance we have two bifur- 



© 2001 RAS, MNRAS 000, §-|l| 



4 Ch. Skokos et al. 



eating families introduced in the system. The number of os- 
cillations along the rotation axis z corresponds to the verti- 
cal resonance at which the family is born. E.g. families xlvl 
and xlv2, which are bifurcated at the vertical 2:1 resonance 
region, have orbits with two oscillations along the z axis. 
This determines only partially their morphology, since the 
bifurcating family can be introduced either in the z or the z 
coordinate of the initial conditions. If we know the number 
of oscillations of a family along each axis and also whether 
it is a bifurcation in z or i, then we know its morphology. 
Families with similar morphology are similar in their corre- 
sponding (x,y), (x,z) and (y,z) projections. In the fiducial 
case, where each vertical resonance is associated with two 
bifurcating families, the families xlv(2n-3) and xlv(2n-2) 
are born at the n:l resonance. 

We note, however, that the first vertical bifurcation is 
not in every model the xlvl family, as in the fiducial case. In 
other models (Skokos, Patsis & Athanassoula, 2002 - paper 
II) it can happen that the first 3D bifurcation of xl is not 
related with the 2:1 vertical resonance, but with a differ- 
ent one. In such a model the first vertical bifurcation of xl 
will have the same name as the family of the fiducial model 
which has similar morphology. Equivalently, it will have the 
same name as the family of the fiducial model which is in- 
troduced in the same n:l resonance and in the same (z or z) 
coordinate. This way we make sure that families with similar 
morphologies share the same name in the various models. In 
addition, if for some reason we have more than one vertical 
bifurcation of xl associated with a vertical resonance, we 
introduce appropriate primes in our nomenclature. E.g. in 
a model with two vertical 4:1 resonances we will have the 
pairs of bifurcating families xlv5, xlv6 and xlv5', xlv6'. By 
keeping the basic name of the family similar for all families 
associated with the same vertical resonance, we underline 
again the dependence of the name on the encountered mor- 
phology. Nevertheless, the basic names are given in the fidu- 
cial model, which thus becomes a reference case for all our 
work. 

We use the same nomenclature not only for the bifur- 
cations of the basic family xl, but in general for the vertical 
bifurcations of every 2D family. Their name consists of the 
name of the parent family, followed by 'vn', where n in- 
dicates its n-th vertical bifurcation. Also the names of the 
bifurcations of the bifurcating families are characterized by 
the addition of '.1', '.2'. . . etc. at the end of the name of the 
3D family, as described above for the corresponding families 
associated with xl. 

We will use the same system in order to name also radi- 
ally bifurcating families. In general a radial bifurcation will 
be named as 'wrn', where 'w' the name of the parent family. 
E.g. the n-th radial bifurcation of family f will be 'frn' (frl, 
fr2...etc). 

Let us now turn to orbits related with the axis of ro- 
tation. The family on the axis of rotation is called 'z-axis' 
family (Martinet & de Zeeuw 1988). Its two first bifurcations 
are introduced at the first S^U transition and the first U-^S 
transition respectively and they are the 'sao' and 'uao' orbits 
of Heisler, Merritt & Schwarzschild (1982). This nomencla- 
ture, however, does not lend itself to extension which can 
include what Poincare (1899) called the 'deuxieme genre' 
families (cf. Polymilis, Servizi & Skokos, 1997) which can 
play an important role in some Ferrers bars (paper II), so 



we will not adopt it here for other families related with the 
z-axis orbits. In practice 'deuxieme genre' orbits are found 
on the stability diagrams as bifurcations of the parent family 
when this family is considered as being of higher multiplic- 
ity, i.e. if its orbits are repeated many times. Thus, the z-axis 
family, when its orbits are repeated twice, is called z2. Bifur- 
cations of the z2 family are called z2.1, z2.2 etc. The same 
rule applies for the bifurcations of z3, i.e. for the bifurca- 
tions of z-axis if this is described three times. We then have 
z3.1, z3.2 and so on. These bifurcating families always come 
in pairs. A further index (s or u) is attached to their names 
and is related with their stability. 

Around the Lagrangian points £4,5 we have the long pe- 
riod banana-like orbits, which form a tree of families, and the 
short period orbits. For the latter we keep the Contopoulos 
& Grosb0l (1989) notation (spo). For the banana- like or- 
bits we use the notation banl, ban2,...bann in the 2D cases. 
Their 2D bifurcations are the families bann.l, bann.2,... etc. 
and their 3D bifurcations the families bannvl, bannv2,... 
etc. 3D banana-like orbits not related with a 2D one are 
named banvn. 

A 2D family found around the unstable Lagrangian 
points 1/2,3 is called l\. 

Throughout the papers we give also the names used by 
other authors for families that have been previously studied. 
However, since our study is more extended, there are several 
families mentioned for the first time. 



3 THE MODEL 

3.1 The 3D potential 

For our calculations we used a 3D potential, which consists 
of a Miyamoto disc, a Plummer bulge and a Ferrers bar. 
Pfenniger and collaborators have made extensive use of this 
potential for orbital calculations (Pfenniger 1984, Pfenniger 
1985a;b, Martinet & Pfenniger 1987, Pfenniger 1987, Pfen- 
niger 1990, Hasan, et al. 1993, Olle & Pfenniger 1998). Our 
work is, in many ways, more extended. We make a much 
more extensive search for periodic families and we further- 
more follow their stability. The latter allows us to find a 
number of 'new' families, which show interesting morpho- 
logical characteristics. Furthermore, we vary the parameters 
of the model so that we are able to make comparisons be- 
tween fast and slow rotating bars as well as between strong 
and weak bars (paper II). Finally, we focus our work more 
on tracing the orbital behaviour that could support observed 
morphological features and less on studying in depth qual- 
itatively the dynamical phenomena that take place in this 
kind of Hamiltonian systems. 

Our general model consists of 3 components. The disc is 
represented by a Miyamoto disc (Miyamoto & Nagai 1975), 
the potential of which reads: 



3>i 



GM D 



vV +y 2 + (A+ VB 2 + z' 2 



(4) 



where Md is the total mass of the disc, A and B are the 
horizontal and vertical scale lengths, and G is the gravita- 
tional constant. The bulge is modeled by a Plummer sphere 
with potential: 
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GM S 



\Jx 2 + y 2 + z 2 + e| 



(5) 



where e 3 is the scale length of the bulge and Ms is its total 
mass. The third component of the potential is a triaxial 
Ferrers bar, whose density p(x) is: 



p(x) = 



105M B 
32Tvabc 



(1 



for m ^ 1 
for m > 1 



(6) 



+ 



a > b > c, 



(7) 



where 

2 2 

2 y . a; 

a 2 o 2 

a, fo, c are the semi-axes and Mb is the mass of the bar 
component. The corresponding potential <l>s and the forces 
are given in Pfenniger (1984)n. They are in a closed form, 
well suited for numerical treatment. For the Miyamoto disc 
we use A=3 and B=l, and for the axes of the Ferrers bar we 
set a : b : c = 6 : 1.5 : 0.6, as in Pfenniger (1984). We note 
that these axial ratios are near the standard values given 
by Kormendy (1982). The masses of the three components 
satisfy G(M D + M s + M B ) = 1. The length unit is taken 
as 1 kpc, the time unit as 1 Myr and the mass unit as 2 x 
10 n M Q . 

In Table |l| we give the parameters of our model. We 
give it the name Al, and it will be one of the models to be 
used in our comparative study in paper II. 



3.2 The 2D Ferrers bar 

The general orbital structure in potentials including a 2D 
Ferrers bar can be found in Athanassoula (1992a). The dy- 
namics are dominated by the presence of the xl family, 
which is in general stable. It is characterized by the pres- 
ence of a narrow instability zone at the 3:1 resonance and 



a gap at the 4:1 region, w hich is generally of type 2 (Con 
topoulos & Grosb0l 1989). The S^U^S transition at the 



3:1 region introduces in the system a couple of simple pe- 
riodic families of orbits, the importance of which remains 
local. Beyond the type 2 gap and above the local maximum 
of the characteristic of xl at the 4:1 resonance (Fig. 2 in 
Contopoulos & Grosb0l 1989) one can find a large number 
of families squeezed close to the zero velocity curve. Finally 
the families x2 and x3 generally exist for a large energy 
range and their characteristics form a single bubble. As it is 
known, x2 is generally stable and x3 unstable. 

In the next sections we describe the orbital behaviour 
in a 3D case where both radial and vertical 2:1 resonances 
exist. We will thus find the differences introduced in the 
morphology and stability of the families of periodic orbits 
by the inclusion of the third dimension. We will also examine 
how the 3D families of periodic orbits support the bar. 



4 THE xl-FAMILY AND ITS BIFURCATIONS 



* We made use of the offer of Olle & Pfenniger (1998) for free ac- 
cess to the electronic version of the potential and forces routines. 




Ej 



Figure 1. xl, x2-x3 and the x2mul2 (Ej , ^-characteristics. The 
curve corresponding to x2mul2 is the projection of its character- 
istic on the (Ej,x) plane. Stable regions are drawn in black and 
unstable ones in light gray. 



4.1 A general description 

Contrary to the 2D models, where a single family, the xl 
family, provides the building blocks for the bar, in 3D mod- 
els we have a tree of families consisting of 2D and 3D families 
related to the planar xl orbits. In Table [] we summarize the 
properties of these families. We list their name, the value 
of the energy at which they are born (Ej), the Ej inter- 
vals where they are stable and we indicate whether they 
are two-dimensional (2D) or three-dimensional (3D). Their 
interconnections and their role will be described in the fol- 
lowing paragraphs. 

There are also several 2D families, which are radial bi- 
furcations of xl and thus part of the xl-tree, but play a less 
important role in the morphology of the models. They are 
described in a separate table (Table ^). The 't' families are 
related with the 3:1 and the 'q' with the 4:1 radial resonance 
region. 

Besides the orbits related to the xl family, we find the 
x2 and x3 families and their 3D relatives as well. They exist 
for the same energy intervals as the families of the xl-tree, 
but their projections on the equatorial plane are elongated 
along the minor axis of the bar. They are described below. 



4.2 Families xl, x2 and x3 

The characteristics of the xl and the x2-x3 families in model 
Al (Fig. [I]) have t he typical geometry o f the characteristics 
of 2D Ferrers bars (Athanassoula 1992a). Due to the vertical 



instabilities, however, xl becomes unstable over several Ej 
intervals, and not only at the radial 3:1 resonance region, as 
in the 2D case. In Fig. [I] and in all characteristic diagrams 
hereafter we draw the unstable regions in light-grey. We ob- 
serve that the decreasing part of the xl curve, below the local 
maximum at the radial 4:1 resonance region (Ej w —0.21), 
is almost everywhere light-grey, indicating that the family 
is unstable there. The curve at about Ej ~ —0.17 turns 
back towards lower energies, remaining after that continu- 
ously unstable. The morphological evolution of the xl orbits 
is the one expected from the 2D case and is given in Fig. [| 
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Table 1. The parameters of our fiducial model Al. G is the gravitational constant, Mo, Mg, Mj are the masses of the disk, the bar and 
the bulge respectively, e s is the scale length of the bulge, Qi, is the pattern speed of the bar, _Ej(r-IILR) and Ej(v-ILR) are the values 
of the Jacobi constant for the radial and vertical 2:1 resonances and R c is the corotation radius. 

GM D GM B GM S e a U b %(r-IILR) F^v-ILR) R c 
0.82 0.1 0.08 0.4 0.054 -0.44 -0.36 6.13 



Table 2. The families of the xl-tree. The successive columns give the name of the families, the value of the energy at which they are 
introduced (Ej), the intervals of Ej at which they are stable and also if they are 2D or 3D. The 'bow'-region is explained in the text, 
while after an energy value indicate that a family continues to be stable, but reaches distances far away from the z = plane. 



family 


E* 

3 


stable intervals in Ej 


2D / 3D 


xl 


-0.495 


-0.495 < Ej < -0.360 


2D 






-0.343 < Ej < -0.293 








-0.278 < Ej < -0.244 








'bow'-region -0.222 < Ej < -0.214 








-0.211 < Ej < -0.205 








-0.192 < Ej < -0.191 








-0.186 < Ej < -0.185 








-0.175 < Ej < -0.173 




xlvl 


-0.360 


-0.360 <Ej < -0.336 


3D 






-0.253 < Ej < -0.147 . . . 




xlv2 


-0.343 


always unstable 


3D 


xlv3 


-0.293 


-0.293 < Ej < -0.221 


3D 


xlv4 


-0.278 


-0.224 < Ej < -0.149 


3D 


xlv5 


-0.213 


-0.213 <Ej < -0.172 


3D 


xlv6 


-0.211 


always unstable 


3D 


xlv7 


-0.205 


-0.205 < Ej < -0.183 


3D 






-0.174 < Ej < -0.170 




xlv8 


-0.192 


always unstable 


3D 


xlv9 


-0.185 


-0.185 <Ej < -0.182 


3D 



The numbers at the upper right corners of the individual 
frames correspond to the Ej value of the orbit. The orbits 
are chosen along the characteristic curve starting from the 
lower values of the Jacobi constant; the orbits in Fig. &L,iJ 
belong to the decreasing branch. Except for the instability 
zones related to the 3:1 resonance all other unstable parts 
of xl appear only in the 3D case. As mentioned in section 
2.1, the families introduced at the instability strips by bifur- 
cation inherit the kind of stability of the parent family, i.e. 
of xl. Thus, the instability gaps on the xl characteristic are 
covered by the stable orbits of the families born after the 
corresponding S— >U transitions. So for almost every energy 
Ej there exists a stable orbit of the xl-tree. As we men- 
tioned in section 2, the 3D bifurcated families are in general 
characterized by four initial conditions (xo, xq,Zo,Zo) so that 
a (Ej,xo) characteristic diagram cannot provide all the es- 
sential information. For this reason we prefer to follow the 
dynamical evolution of the orbits using stability diagrams. 
These diagrams frequently become complicated, but they 
have the big advantage of giving in a straightforward way 
the interconnections of the various families, thus becoming 
a very useful tool in the hunting of periodic orbits. 

The evolution of the stability indices b\ and &2 for xl are 
given in Figs. ^, ^, and ^for successive energy intervals. The 
arrows denote bifurcated families at the bifurcating points 
and show the direction of the stability index associated with 



bl,2 
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jy \xivi 
- — ^ V) 
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/ VL^--t\ x1 ti 
v / \/ 




\ VA \ 

xlv2 xlv4 \ 



-0.5 -0.45 -0.4 -0.35 -0.3 -0.25 



Ej 

Figure 3. First part of the xl stability diagram. Arrows denote 
the bifurcation of families and the direction in which these evolve. 

the S— >U or U^S transition. We observe that the variation 
of the index which in Fig. ^| has the larger values for Ej < 
—0.38 brings in the system the 3D families xlvl, xlv2 etc., 
while the variation of the other index brings in the families 
associated with the radial instabilities. The latter remain on 
the equatorial plane. The variation of their stability indices 
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Table 3. Radial bifurcations of the xl family, 't' families are related with the 3:1 and 'q' families with the 4:1 radial resonance region. 
Columns are as in Table ^. Ej values are given in general with three digits, except from the cases where narrow Ej ranges of existence 
need more accuracy. We note that the t3 family exists also for lower energies than its E* (see section 4.4). 



family 


E*. 
j 


stable intervals in Ej 


2D / 3D 


tl 


-0.244 


-0.244 <Ej < -0.218 


2D 


t2 


-0.214 


-0.214 < Ej < -0.209, 


2D 






-0.204 < Ej < -0.203 




t3 


-0.205 


-0.2065 < Ej < -0.2005 


2D 


ql 


-0.191 


always unstable 


2D 


q2 


-0.1857 


-0.1860 < Ej < -0.1857 


2D 


q3 


-0.183 


-0.1818 < Ej < -0.1808 


2D 
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Figure 4. Second part of the xl stability diagram. The evolution 
of the stability curves in the xl 'bow' is indicated with numbers 
from 1 to 8 and arrows. In the diagram we indicate also the xl 
bifurcations to the right of the 'bow'. They are the families t2, 
xlv5, xlv6 and xlv7. 



will in turn bring new families due to vertical and radial 
instabilities. 

The feature depicted in Fig. ^ is typical of the stability 
diagrams of many of our models. We call this kind of evolu- 
tion of the stability indices a 'bow'. The bi and 62 curves do 



not break anywhere, but they evolve in a continuous, rather 
complicated way, changing direction twice. This 'bow' area 
corresponds to the bend, or elbow, in the characteristic at 
about E 3 « -0.227 (Fig. [j]), and the complicated evolution 
of the stability indices happens as we move towards lower 
Ej values along the characteristic curve of xl at this area. In 
Fig. [j] one can follow the evolution of b\ and i>2 by following 
the evolution of both the numbers and the nearby arrows. 
The lowest value of the stability index at '5', not included 
in the Figure (indicated only with a dashed arrow outside 
of figure frame), is « —55. 

A significant change in the way the 3D bifurcations of 
xl are introduced in the system happens at the instability 
zone found just beyond the local maximum of the (Ej,x) 
characteristic close to the radial 4:1 resonance. As we see in 
Fig. |3j ^ and [|, the 3D families are bifurcated at S— >U and 
U — >S transitions, where the corresponding stability index 
intersects the b = — 2 axis. Moving on the characteristic 
towards corotation, before reaching the decreasing branch, a 
bifurcating family at an S-^U transition is a stable 3D family 
with initial conditions (x,z,x,z) = (a, b, 0,0), where a, b 6 
R and a,b 7^ 0. On the other hand, the family bifurcated 
at the U— >S transition, is (initially) simple unstable and 
has initial conditions (x,z,x,z) = (c, 0,0, d), with c, d £ R 
and c, d ^ 0. This means that the family introduced in the 
system as stable is a bifurcation at 2, and the simple unstable 
family a bifurcation at i. For the set of families associated 
with the vertical 5:1 resonance, on the decreasing branch 
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Ej 



Figure 5. Stability diagram for xl and bifurcating families cor- 
responding to the decreasing part of the xl characteristic. The 
evolution of the stability curves is indicated with numbers from 
1 to 11, and thick arrows that point to the direction of the evolu- 
tion. The bifurcating families and their direction of evolution are 
denoted with thin arrows. 



of the characteristic, this sense of bifurcation is reversed. 
Namely we have the bifurcation in z at the U— >S transition 
(xlv8) and the bifurcation in i at S— >U (xlv7). 

In Fig. ^|we plot the last part of the stability diagram of 
the xl family, corresponding to energies higher than —0.2. 
As can be seen from the characteristic diagram of Fig. this 
includes most of the decreasing part of the characteristic, the 
bend at Ej ss —0.173 and the part that goes towards lower 
energies. This part (roughly for —0.22 < Ej < —0.173), has 
negative x values starting soon after the bend. Heavy ar- 
rows and numbers in increasing order on and next to the 
stability curves in Fig. [| indicate the evolution of the in- 
dices as we move along this part of the characteristic. As 
we can see most parts are unstable, the short stable parts 
being drawn with heavy lines. After the turning point, at 
Ej w —0.173, the upper curve, moving now towards lower 



Ej values, is stable until Ej 



-0.181, then has a part with 



values smaller than —2 and then reenters the stability re- 



gion for Ej 



-0.197. The lower stability curve, however, 



reaches absolutely large negative values. Thus, the family 
is always unstable in the parts where x < 0. It is easy to 
understand how the negative x values are introduced by fol- 
lowing the evolution of the xl orbit morphology as we move 
along the characteristic (Fig. |^). As one moves along the 
decreasing part of the characteristic (Fig. ^i — >j), the four 
apocentra of the orbits develop loops, whose size increases 
strongly as the energy increases. Already for the orbit in 
Fig. fcjj the loops have become so large, that the sides of the 
orbit along the bar major axis nearly touch. As we continue 
along the characteristic they will touch and then cross, so 
that x becomes negative. 

Let us now present the evolution of the x2-x3 loop in 
the 3D case. As seen in Fig. the situation with the x2-x3 
characteristic is exactly like in 2D. The stability indices form 
also loops, as the b\ and 62 indices of x2 and x3 join each 
other in pairs (Fig. []). In 2D models, the families xl and x2 
are in general the only simple periodic stable families at the 
x2-x3 area. This is not necessarily the case in 3D models. 




E, 



-0.36, while close to E-j 



becomes again stable at about Ej 



xlvl as stable up to Ej 



-0.147. 



-1 -0-5 05 1 -1 -05 0-5 1 -1 -0-5 05 1 

Figure 7. The 3D x2-like orbit x2mul2. This is a family of mul- 
tiplicity 2. 



E.g. in this model, as we can see in Fig. gb, the 3D fam- 
ily xlvl has been bifurcated as stable just before the point 



-0.29 the family xlv3 is 



introduced in the system. So the situation at the x2-x3 area 
is more complicated, since we have there four simple peri- 
odic stable families. Since the x2-x3 stability indices form a 
bubble they have no further intersections with the b — —2 
axis and there are no further bifurcations of other simple 
periodic x2-like families. Both families, however, have tan- 
gencies with the b = 2 axis. At these points, as mentioned in 
the introduction, families of the same kind of stability, but 
with double multiplicity, will be bifurcated. The one bifur- 
cated from the stable family x2 is interesting. If we put its x 
initial values on the characteristic diagram (Fig. ^), we ob- 
tain the extra branch emerging from the x2-x3 loop, pointed 
with the curved arrow and characterized as 'x2mul2'. The 
energy range over which it is stable is indicated with a dou- 
ble arrow above the 6 = 2 axis in Fig. Its morphology 
is given in Fig. The (x,y) projection is typical of an x2 
orbit, the (x, z) one is a fish-like figure reflecting the double 
multiplicity of the family, while the (y, z) projection offers a 
shape that could produce a tiny boxy structure in the cen- 
tral region of the bar (note the scale on the axes). The (x, z) 
projection can also offer a boxy structure, if one considers 
together with every orbit its symmetric with respect to the 
z axis. This, however, is elongated along the minor axis of 
the bar as will be discussed in paper III. The morphology 
of this family shows that the model clearly can support in 
its face-on projection the presence of stellar rings in the x2- 
x3 area. This, however, is a thick ring structure extending 
outside the equatorial plane. 



4.3 The main 3D families 

The xl S^U transition at about Ej ^ -0.36 (Fig. |b), 
generates the 3D family of periodic orbits xlvl. This family 
is related with the presence of the vertical 2:1 resonance. 
It has a stable part close to the bifurcating point, then it 
has a complex unstable part after an S— >A transition, and 



-0.253. We have found 



The morphological evolution of the family xlvl is given 
in Fig. || This family corresponds to the z2 family of Hasan 
et al. (1993) and its orbits have been associated with the 
appearance of the peanut shaped bulges by Combes, Deb- 
basch, Friedli et al. (1981). Indeed, due to the symmetry of 
the potential with respect to the equatorial plane, one can 
find all 3D families in pairs. Thus for xlvl we will have the 
smile (^) and frown (<->) types of the (y, z) edge-on projec- 
tions coexisting at a given energy, and the same holds for 
the (x, z) projection. The (x, y) projections of the 3D or- 
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Figure 6. Stability diagram for xl, x2 and x3 orbits. In order to follow the interconnections of the various families, this diagram is given 
in two parts. In (a) the stability indices of the families x2 and x3 are emphasized and in (b) those of xl and of the bifurcating families. 
A horizontal segment with double arrows in the upper part of the diagrams indicates the range of stability of the family x2mul2. The 
horizontal segment with double arrows, drawn black in the lower panel and indicated with A, denotes the complex unstable part of xlvl. 
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Figure 8. Three stable orbits of the xlvl family. Note that the 
upper panels have a different scale than the middle and lower 
ones. Corotation in model Al is at 6.13. The energies from top 



to bottom are: E; 



-0.35, -0.25 and -0.20 respectively. 



bits follow in general the morphology of the corresponding 
xl orbit of the same energy. As we said in the introduc- 
tion, the importance of a family of stable periodic orbits is 
limited as the individual orbits grow in \z\. The xlvl orbit 
for Ej — —0.2 in Fig. exceeds both in its (x,z) and its 
(y, z) projections the height of 2 kpc and this means that it 
cannot contribute significantly to the density of the galactic 
disc. Its spatial extent on the other hand indicates that this 
orbit could be used to populate the bulge area. 



The U— >S transition at Ej 



-0.343 generates the fam- 



ily xlv2 which we followed until Ej « —0.173. It remains 
totally unstable and ends after a U — > D — > A sequence. It 
thus doesn't play any important role in the dynamics of the 
system. 

Family xlv3 (Fig. ^) is stable and its orbits keep low \z\ 
values roughly in the interval —0.293 < Ej < —0.221. It then 
ends with a S— >A transition. This family is similar to the 
zl family of Hasan et al. (1993). We note that both xlvl 
and xlv3 provide useful orbits in the system before their 
S^A transition. This behaviour is also seen in the 3D thick 
spiral model in Patsis & Grosb0l (1996). Complex instabil- 
ity helps introducing abrupt drops in the density of given 
features of a model (in our case the peanut), since it stops 
abruptly the existence of the family responsible for their ap- 
pearance without bringing new stable families in the system. 
On the other hand, in cases where a stable family donates 
its stability to a bifurcation we have a smooth morphologi- 
cal evolution, which can give smooth density profiles in the 
galaxies. Both xlvl and xlv3 do not have any intersections 
or tangencies with the —2 axis and for this reason they do 
not bifurcate other families with the same multiplicity. 

The next bifurcated family is xlv4. This is bifurcated 
from xl after a U— >S transition. We would thus have ex- 
pected it to be unimportant, since its parent family, xl, is 



Figure 9. Three stable orbits of the xlv3 family. Their (x, z) and 
(y, z) projections have always low \z\ values. The energies from 



top to bottom are: Ej 



-0.28, -0.26 and -0.22 respectively. 



unstable at the bifurcating point. This is the typical be- 
haviour in such cases and we have seen it already happening 



for xlv2. xlv4 is introduced at about Ej 



-0.278. One 



of the two stability indices, let us call it &i, remains in the 
interval —2 < bi < 2, while the other, &2, goes to negative 
values smaller than —2. For larger Ej values, however, &2 
increases and for about Ej ~ —0.224, both indices come 
in the stability zone, i.e. we have —2 < 61,2 < 2. The de- 
tailed description of this complicated evolution is beyond the 
scope of the present paper and does not add anything to the 
important information that family xlv4 brings stable repre- 
sentatives in the system for Ej > —0.224. The xlv4 family 
remains stable up to Ej ~ —0.149 where it becomes simple 
unstable. Its stability indices fold and the family continues 
existing towards smaller energies. The morphological evolu- 
tion of xlv4 can be seen in Fig. In Fig. |Io| a we give the 
three projections of an unstable orbit close to the bifurcating 
point from which the family emanates, while in Fig. p"o| b and 
Fig. [[fjfc we give two stable orbits, for energies Ej > —0.224. 
The last one is for Ej — —0.206 and we see that already 
the orbit reaches \z\ values close to 2 kpc away from the 
equatorial plane. For each orbit of this family there is also a 
symmetric one with respect to the equatorial plane. If only 
one of the two is populated, this would give rise to an asym- 
metric warp-like shape. Populating them both restitutes of 
course symmetry. The stable orbits of this family enhance 
the bar, but they deviate substantially from the equatorial 
plane. 



4.4 Families at the 3:1 radial resonance 

As we already saw, one of the two stability indices bifurcated 
the 3D families xlvl, xlv2, xlv3 and xlv4, by its intersec- 
tions with the b — —2 stability axis. The intersections of 
the second stability index with this stability axis introduce 
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Figure 10. Orbits of the family xlv4. Panels (a) show an unstable 
orbit close to the bifurcating point at Ej = —0.278. Panels (b) 
and (c) show stable orbits for Ej 
respectively. 
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Figure 11. Stability diagram of the tl family, the first radial 
bifurcation of xl at the 3:1 resonance. The stability indices of xl 
are given as well, drawn with light grey lines. 



in the system planar 2D orbits. The first family is bifur- 
cated after a S^U transition at Ej w —0.244, i.e. in the 3:1 
resonance region. We call it tl and it is stable (Fig. [H]). It 
bridges exactly the instability zone of the xl in the S^U^S 
transition, i.e. its s tability indices toge ther with those of the 
xl, form a bubble ( Contopoulos 198q ). tl exists for approx- 
imately -0.244 < Ej < 



-0.218 and at Ej « -0.218 can be 



considered as an inverse bifurcation^] of xl. At Ej ~ —0.214, 



t Inverse bifurcation is a non-linear phenomenon encountered in 
Hamiltonian systems, according to which the bifurcated family, 
instead of evolving towards the same direction as the parent fam- 
ily, changes direction. It thus extends for the same energies as the 




Figure 12. Stable orbits of the three 2-dimensional families at 
the 3:1 area. Note that the loops of t3 are asymmetric. In the left 
panel we plot a tl orbit together with its symmetric with respect 
to the bar minor axis. 



just beyond the 'bow' area, the same stability index has an- 
other intersection with the b = — 2 axis and xl bifurcates 
another 2D family, t2. Several 2D and 3D 3:1 type fami- 
lies, related to each other and with xl, are introduced in 
the interval —0.214 < Ej < —0.20. Let us briefly mention 
that, besides tl (in Fig. |Tl| ) and t2, we found a third 2D 
3:1 family, t3, which is stable for -0.2065 < Ej < -0.2005, 
although it is introduced in the system as simple unstable 
for Ej w —0.205. The morphology of the three 2D families 
tl, t2 and t3 is given in Fig. [l2] and their stable energy in- 
tervals in Table pi For the lower energies, the tl orbits have 
only one loop, which is located on the y axis, as the exam- 
ple shown in the left panel of Fig. [l2| For higher energies 
they develop two more loops, symmetric with respect to the 
y axis, and roughly equal in size to the first one. Since the 
orbits of family tl are symmetric with respect to the y axis, 
for every orbit we should have also its symmetric with re- 
spect to the x axis. Combining the two, as in the left panel 
of Fig. |l2| we obtain a shape that is elongated along the bar 
major axis and resembles the morphology of the xl orbits 
with loops, at least for the energies where the orbits have 
only one loop. The extent of such orbits along the y axis 
reaches up to 4 kpc, i.e. two thirds of the way to corotation. 

t2 brings in the system three-dimensional families of 
periodic orbits with stable representatives. It bifurcates the 
family t2vl at Ej ~ —0.209, which in turn bifurcates t2vl.l 
at Ej w —0.205. The t2vl family provides stable orbits to 
the system for -0.209 < Ej < -0.207 and the t2vl.l family 
for -0.205 < Ej < -0.203. Triangular-like t2-type orbits 
have characteristic peaks at the sides of the bar, like the 
peak of the orbits at x w —4 in the (x,y) projections of 
Fig. They are near but not always on the minor axis of 
the bar and their presence can lead to local enhancements of 
the density at the area between the bar and the L,4 ; 5 points. 
For any energy in the interval —0.214 < Ej < —0.20 there 
are almost always stable 3:l-type orbits of one or the other 
family. Together with tl, they affect the dynamics of the bar 
in this region. We note that the 3:1 orbits bifurcated from 
xl are very common in all barred potentials and have both 
in 2D and 3D dynamically only local importance. Orbits of 
type t2vl and t2vl.l have been found even in the early N- 
body simulations of 3D bars (Figure 5 in Miller & Smith 
1979). The loops of t3 on either side of the major axis of 
the bar are not of equal size. As can be seen by careful 
inspection of the t3 orbit in Fig. |l^, the loop at the right side 



parent family before the transition and h as the kind of stab ility 
of the parent family after the transition. (Oontopoulos 198q). 
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Figure 13. Stable orbits of two 3-dimcnsional families at the 3:1 
area. Their names are given on the top left of each sets of panels. 




Figure 14. A typical stable xlv5 orbit. 

of the major axis is slightly bigger than the one to the left. 
Thus, morphologically, t3 is a kind of asymmetric tl, since 
for larger energies tl develops loops which are symmetric 
with respect to the major axis, besides the one along the 
major axis. 

4.5 The last part of the xl-tree 

There are two more 3D bifurcations of xl close to the lo- 



cal maximum of the characteristic at 



-0.205. It is 



xlv5 (bifurcated at Ej 



-0.213 and being stable un- 



til Ej « -0.172) and xlv7 (bifurcated at Ej w -0.205, 
just beyond the peak of the characteristic). The family 
xlv7 and its bifurcation xlv7.1 provide stable orbits for 
-0.205 < Ej < -0.18 and -0.175 < Ej < -0.17. Neverthe- 
less, the part of this family that contributes to the density of 
the bar is limited by the fast increase of \z\ with the energy. 
Fig. [H| and ^ show the morphology of these families. 

In the same region we encounter two more 3D bi- 
furcations of xl, namely the families xlv6, introduced at 
Ej » -0.211 (Fig. [J), and xlv8 introduced at Ej » -0.1925 
(Fig. ^). Both are born after an U— >S transition of xl and 
remain always unstable. We note that the representatives 
of xlv5, xlv6, xlv7 and xlv8 families are morphologically 
similar to those of the Bz2, Bi2, B23 and BZ3 families of 
Pfenniger (1984) respectively. 

As we have seen, xl is mostly unstable in the decreasing 
branch beyond the local maximum at the radial 4:1 gap and 
the morphology of the orbits at this branch is in general 
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Figure 15. (a) A xlv7 orbit and (b) a xlv7.1 one. Both are 
stable. 
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Figure 16. Morphology of stable orbits of xl's bifurcations at 
largest energy values. Panels (a) and (b) show members of families 
q2 and q3 respectively. Panels (c) to (e) show the three views of 
an orbit of the xlv9 family. 



rectangular-like with loops in the corners. There are several 
families bifurcating from this branch and their orbits have, 
as already noted for other families, a morphology similar to 
that of xl in the region. The 2D families q2 and q3 provide 
stable asymmetric orbits, two examples of which are given 
in Fig. hq. We also have one 3D bifurcating family, xlv9, a 
member of which is shown in Fig. This family also has an 
asymmetric stable bifurcation for a short energy interval. No 
stable members of these families can be found outside the 
interval -0.186 < Ej < -0.1808. To this we should add the 
small intervals of stability provided by xl itself (cf. Fig. [j] 
and Fig. |). 

Finally, for the sake of completeness we give in Fig. ^ 
the morphology of the three 3D families, members of the xl- 
tree, which remain always unstable although they exist for 
large energy intervals. As we have seen in the corresponding 
paragraphs they are the families xlv2, xlv6 and xlv8. 



5 FURTHER FAMILIES 
5.1 Orbits around L4 and L5 

Another important 'forest' of families is the group of the 
banana like orbits. Here we find the usual planar long and 
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Figure 17. Orbits of the unstable families xlv2 (a), xlv6 (b), 
and xlv8 (c) at Ej = -0.2612, -0.178 and -0.184 respectively. 



sho rt p eriod banana-like orbits (Contopoulos & Grosb0l 
198S |). j rhe long period orbits are coming in the system in 
a large variety of families all of which have stable parts for 
-0.1984 < Ej < -0.1944. The stability indices of these 
orbits exhibit a complicated behaviour having several tan- 
gencies and intersections with the b — 2 and b = —2 axes. 
This brings many families in the system by bifurcation. The 
family found for lowest Ej values is banl (Fig. IS a) which 
is born at Ej « -0.1984, followed by ban2 (Fig. lib,c) that 
appears at a slightly greater energy value - which in turn 
bifurcates ban2.1 (Fig. |l|d) at Ej « -0.1972 - and ban3 
(Fig. |l|e,f) introduced in the system at Ej ^ -0.1983. The 
most important of the planar orbits with stable parts, is 
ban4 (Fig. [jj|h), because it is stable over the largest en- 
ergy interval (-0.1982 < Ej < -0.1955). It exists for Ej > 
—0.1982 and it is not bifurcated at this point from any of 
the families existing for lower energies (banl, ban2, ban2.1, 
ban3, ban3.1). From ban4 bifurcates the 2D family ban4.1 
(Fig. |li|g). The stability indices of ban4 have a compli cated 



behaviour which is typical of a collision of bifurcations ( Con 

topoul c|s 198ej )Fl Approaching Ej = -0.1955037765, the 
ban4 orbits shrink to L4 (or L5), and beyond this point the 
short period orbits (spo) grow in size and take their bean- 
like shape (Fig. |li| h and i respectively) . 



t Collisions of bifurcations happen when both f>i and 62 are ex- 
actly equal to —2 or 2 for a particular set of the control parame- 
ters. In order to observe a collision we need to vary continuously 
a control parameter of our model (i.e. to consider successive indi- 
vidual models), and for all these cases to follow the evolution of 
the stability indices as a function of Ej. This practically means 
that we vary two parameters. If it happens that 61 = hi = —2 (or 
2) for a critical set of the control parameters, then we will observe 
a change in the interconnections between parent and bifurcating 
families, before and after the collision. This may also change the 
general behaviour of the dynamical system. 
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Figure 18. Stable 2D banana-like orbits. Ej is given in the lower 
left corner of panels (a) to (g). Panels (h) and (i) include orbits of 
many Ej values. Arrows indicate morphological evolution of the 
same or related families. 



We have also found three 3D families of periodic or- 
bits with stable parts. ban3vl (Fig. [l9^,), a bifurcation of 
ban3 at Ej ~ —0.1982, is initially marginally stable, having 
one of its two stability indices almost equal to —2, but for 
Ej > —0.1962 the index become clearly larger than —2. At 
Ej « —0.1947 the two indices join each other and we have 
a S^A transition. ban4vl (Fig. [l9|b), a bifurcation of ban4 
at Ej w —0.1976, is almost everywhere marginally stable in 
the interval -0.1976 < Ej < -0.1944. For Ej > -0.1944 
it is always complex unstable. ban3vl and ban4vl extend 
to very large Ej values, but as complex unstable. Since we 
have S— >A transitions there are no bifurcating families and 
this is the mechanism that terminates the trapping of mate- 
rial around banana-like orbits in our 3D bars. Finally banvl 

-0.1957 as 



(Fig 



is introduced in the system at Ej 



stable and remains stable up to Ej « —0.1954. This family 
is not obviously related to any other banana- like orbit. Since 
it is a 3D family we name it banvl. 



5.2 Orbits around Li and L2 

The Li , L/2 Lagrangian points ar e known to be always un- 
stable (Binney & Tremaine 1987). Around them we find a 



family of planar periodic orbits we call l\. It appears for 
Ej values larger than the one corresponding to Li , the mor- 
phology of its orbits resembles that of the spo orbits rotated 
by tt/2, and their periods are of the order of the epicyclic 
period. Close to the Li energy and for Ej < —0.168 these 
orbits are unstable. For Ej > —0.168, however, l\ has both 
stability indices between —2 and 2 and the family becomes 
stable. Orbits of this family can be found only by starting 
with initial conditions on the major axis of the bar. For this 
reason they had not been previously found, since in previous 
studies searches for periodic orbits started only with initial 
conditions on the y = axis. In Fig. |20j, we plot a few stable 
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Figure 19. Stable 3D banana-like orbits. All three extend to 
large Ej values but the two most important (ban3vl, ban4vl) 
are at these large Ej values complex unstable. The numbers at 
the bottom of the (x, z) projections give the Ej of each orbit. 




Figure 20. Stable orbits of the l\ family and their symmetries 
with respect to the x axis. The innermost orbit, just next to the 
arrow, corresponds to Ej m -0.168, just after the U— >S transition. 



orbits of l\ and their symmetries with respect to the x axis 
for Ej > —0.168. These stable orbits do not support the bar 
since they are elongated parallel to the minor axis. Neverthe- 
less, they are of physical interest since they support motion 
parallel to the minor axis, contribute to the exchange of ma- 
terial between regions inside and outside corotation and are 
able to influence the dynamics in the region between bar and 
spirals in barred spiral galaxies. The streaming at the apoc- 
entra of the l\ orbits could support arc-like features beyond 
the end of the bar. 

For larger energies the l\ orbits can be observed shifted 
towards the x axis (minor axis of the bar), at about Ej ~ 
—0.12 they cross the x axis and after a short unstable zone 
they fall on the retrograde family x4 as stable. 



5.3 Orbits outside corotation 



Beyond corotation we find the usual planar families (Con 
topoulos fc Grosb0l 1989). Most of their members display 
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Figure 21. A 3D stable orbit beyond corotation. 




Figure 22. Examples of stable 2D orbits beyond corotation. They 
support motion parallel to the minor axis of the bar at the coro- 
tation region and favour exchange of material between regions far 
from the center and regions near to it. 



loops. We also find several 3D families with stable parts. As 
an example we give the family depicted in Fig El| which is a 
bifurcation of the planar family called xl(l) by Contopoulos 
& Grosb0l (1989). The vertical extent of the 3D orbits we 
found beyond corotation is in general small. 

Let us also mention some 2D families, orbits of which 
are given in Fig. They have been calculated starting with 
initial conditions on the major axis of the bar as in the case 
of the i\ family and have thus not been described in previ- 
ous papers. All of them have large stable parts. These orbits 
have two interesting properties. First, some of them could 
support motion close to the end of the bar parallel to its 
minor axis, at radii shorter than the corotation radius. Sec- 
ond, they could efficiently transport material from the outer 
parts of the disc, e.g. from a distance close to 20 kpc from 
the center, to the central regions of the bar (e.g. Fig. |22|a) . 
This is particularly true for orbits as those shown in Fig~|22fi 
and c. 



6 CONCLUSIONS 

In this paper we have made an extensive study of both the 
2D and 3D periodic orbits in a fiducial model representative 
of a barred galaxy. We report on the stability and morphol- 
ogy of the main families. Our main conclusions are: 

(i) So far the xl orbits were considered the backbone of 
bars. This, however, can only be the case for 2D bars, since 
the xl can only populate the z=0 plane. For 3D bars the 
backbone is the xl, together with the tree of its 3D bifur- 
cating families. Trapping around these families will deter- 
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mine the thickness and the vertical shape of galaxies in and 
around the bar region. Major building blocks for the 3D bars 
can be supplied also by families initially introduced as un- 
stable. Thus the family xlv4, introduced in the system after 
a U^S transition, is a basic family, giving stable represen- 
tatives for large energy intervals in the system. 

(ii) The (x, y) projections of the 3D families of the xl-tree 
retain in general a morphological similarity with their parent 
family at the same energy. This has important implications 
for the morphology of a galaxy since it introduces building 
blocks which have similar morphology as the xl orbits, but 
have a considerable vertical extensions. Especially at the re- 
gions close to the bifurcating points the (x,y) morphology of 
a xlvn family is not only geometrically similar, but actually 
very close to the morphology of the corresponding xl orbit. 

(iii) The way the 3D families of the xl-tree are introduced 
in the system at an instability strip determines the impor- 
tance of the bifurcations in z or z. Particularly in the present 
model, all 3D families of the xl-tree at the increasing part of 
the characteristic which are bifurcated in z are introduced 
in the system as stable. On the other hand the stable family 
associated with the 5:1 vertical resonance (xlv7), bifurcated 
at the decreasing part of the xl characteristic, beyond its lo- 
cal maximum, is bifurcated in z. Whether the stable family 
of the last S— >U— >S transition is the bifurcation in z or z 
determines in a large degree the model's morphology at its 
outer parts. 

(iv) The radial 3:1 resonance region provides in the sys- 
tem several 2D and 3D stable families. Their role, however, 
is locally confined, as in 2D models. 

(v) 3D orbits elongated along the minor axis of the bar 
can be given by bifurcations of the planar x2 family. 

(vi) We have found several families of 3D banana-like or- 
bits around L^s. Their extent is always restricted by a S— >A 
transition. 

(vii) Stable families found beyond corotation circulate 
material between the outer parts of the system and regions 
as far inwards as 1 kpc. This contributes to the mixing of 
the elements in a disc galaxy. 

The families of periodic orbits we described up to now 
are indeed the basic families of a 3D Ferrers bar. As we ex- 
plore the parameter space, however, their properties change, 
while new important families may appear and play a crucial 
role. A notable example is z3.1s, a family related to the z- 
axis orbits along the rotational axis, which will be described 
in paper II. However, these are rather particular cases and 
are not encountered in every model. 
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